Stochastically sustained population oscillations in high-/? nanolasers 
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Non-linear dynamical systems involving small populations of individuals may sustain oscillations 
in the population densities arising from the discrete changes in population numbers due to random 
events. By applying these ideas to nanolasers operating with small numbers of emitting dipoles 
and photons at threshold, we show that such lasers should display photon and dipole population 
cycles above threshold, which should be observable as a periodic modulation in the second-order 
correlation function of the nanolaser output. Such a modulation was recently reported in a single- 
mode vertical-cavity surface-emitting semiconductor laser. 

PACS numbers: 42.55.Sa, 42. 60. Mi, 42.50.Ar 



I. INTRODUCTION 



Nanolasers, such as semiconductor microcavity lasers 
[iHl] or plasmonic lasers 041], display cavity volumes of 
the order of a fraction of a cubic wavelength and thus 
contains few emitters, while their lasing threshold is at- 
tained with very few photons in the cavity, of the order 
of /3~^^^ 0] at any given time, where /3 is the fraction 
of spontaneous emission funnelled in the "useful" mode. 
In a nanocavity, /3 may be close to 1, so that the num- 
ber of photons at threshold is very small, typically 1 to 
10. For such small numbers, fluctuations are large (of 
the order of the average values), occur only in integral 
multiples of a base value (corresponding to a change of a 
single unit), are asymmetric (as excursions into negative 
values are impossible) and cause the nanolaser output to 
deviate from that predicted by the traditional laser rate 
equations [8|, |9|. The reason for this deviation is that 
these real-number differential equations have been devel- 
oped for conventional macroscopic lasers, in which the 
discrete nature of the number of emitters and photons 
may be ignored (termed "the thermodynamic limit" 0) 
and are therefore inadequate for describing the operation 
of nanolasers. 

Theoretical methods based on full quantum electrody- 
namic treatments have been developed [lol | to account 
for phenomena observed in single-atom or few-atom mi- 
crolasers and micromasers, which have been studied ex- 
perimentally since the 1980s pT| - [l3| . However, in these 
systems the emitting atoms are relatively well isolated 
from their environment, have sharp optical spectra, and 
generally couple strongly to the cavity mode. Semicon- 
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ductor (or plasmonic) nanolasers, on the other hand, 
operate at the opposite limit: The emitters (generally 
semiconductor quantum dots) are embedded in the cav- 
ity material and thus their radiating dipoles interact ex- 
tensively with their environment, undergoing very fast 
dephasing processes. At the same time, under the strong 
pumping conditions required for lasing, the emitters are 
usually highly-excited, containing several electron-hole 
pairs (dipoles) and producing very broad emission spec- 
tra [14:]. These conditions are contrary to the assump- 
tions of the few-atom microlaser theories, as they pre- 
clude strong coupling between the emitters and the cavity 
mode, while the incoherent fluctuation processes cause 
decay of the atomic polarization and rapidly wash out 
any phase- related phenomena. In this limit, the full 
quantum equations reduce to the traditional laser equa- 
tions that involve only photon and dipole numbers 0- 
However, these are continuous-variable differential equa- 
tions that describe the evolution of the mean photon and 
dipole populations, but ignore their discrete nature: In 
each individual realization of the lasing process only in- 
teger values of photon and dipole numbers may be ob- 
tained, a feature that may be ignored in large lasers but 
becomes important in nanolasers where these numbers 
are small. As described below, the traditional laser equa- 
tions can be adapted to the discrete nature of the pho- 
ton and dipole populations through discrete simulation 
methods in which small increments in time are made, 
and within each increment a discrete but random num- 
ber of constituents is chosen, obeying the statistics of the 
physical process being modelled. 

This procedure is quite similar to that used in mod- 
elling small predator-prey or epidemiological ecosystems 
[3 EB], where it has been shown that "demographic 
stochasticity" , that is population fluctuations due to ran- 
dom birth, death, predation or contagion events, pro- 
duces cyclical population variations. These cycles are 
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due to the discrete nature of the individuals in the pop- 
ulation and cannot be accounted for through traditional 
deterministic models for the populations, based on real- 
number continuous differential equations, such as the 
Lotka-Voltcrra equations, whose domain of validity cor- 
responds to the cases in which discretization may be ig- 
nored, such as in infinite ecosystems, or when extensive 
ensemble-averaging is possible. 

The traditional method of dealing with fluctua- 
tions in lasers consists of inserting randomly-fluctuating 
Langevin forces in the laser rate equations and lineariz- 
ing about the steady-state photon and dipole populations 
[TSi] . However, this method is not suitable when small 
populations are involved, since in that case the magni- 
tude of the fluctuations is of the same order as the mean 
values (and thus linearization is not an adequate approx- 
imation) while, at the same time, the small-population 
statistics deviate signiflcantly from the Gaussian statis- 
tics of the Langevin forces. Thus, while the traditional 
treatment may be adjusted so that the variance of the 
Langevin forces reproduces the observed variance of the 
fluctuations, it cannot account adequately for the higher- 
order rnoments of their distribution. In a recent publi- 
cation Q we measured the higher-order intensity fluc- 
tuation statistics of a pulsed nanolaser (up to the fourth 
order) and, through the implementation of a discrete sim- 
ulation method, we showed that the discrete nature of the 
photon number at threshold induces a jitter in the pulse 
timing, which accounts for the observed photon statistics 
to all orders. 

In this paper, we examine a continuously pumped 
nanolaser and show that the "granularity" of the pho- 
ton and dipole populations and the "discretization noise" 
it introduces lead to large self-sustained oscillations in 
these populations, often reaching extinction periodically 
for the case of photons, even above threshold. While this 
behavior is similar to that of small predator-prey ecosys- 
tems, it can also be understood within the traditional 
laser theory as arising from the laser's relaxation oscil- 
lations driven by the discrete jumps associated with the 
spontaneous and stimulated emission events. The pho- 
ton population cycles in nanolasers should be observable 
as a periodic modulation of g'^^^r), the intensity auto- 
correlation function of the nanolaser. 

Such a periodic modulation has recently been reported 
by Takemura et al. l?'] in experiments involving a 
single-mode vertical-cavity surface-emitting laser (VC- 
SEL). While VCSELs are not nanolasers, the sample ex- 
amined displayed (3 w 10^**, corresponding to /3~^/^ sa 
100 photons in the cavity at threshold, so that a periodic 
modulation of g^^^r), albeit of small amplitude, could be 
measured in the vicinity of the threshold. The discretiza- 
tion noise model accounts well for this modulation, thus 
validating its use in solid-state nanolasers. 



II. DISCRETIZATION NOISE 

We consider an optical cavity with one only "useful" 
(or reference) mode, while all other modes are considered 
as "leaks". The cavity contains an ensemble of N ex- 
cited emitters (dipoles) each emitting photons at a spon- 
taneous rate 7||. We assume that the emitters interact 
extensively with their material environment, so that they 
undergo very rapid dephasing processes (as is the case, 
for example, in solid-state room-temperature lasers) in a 
time short compared with the characteristic time of the 
coupling with the electromagnetic field. 

Under these assumptions, the quantum mechanical 
Maxwell-Bloch equations (or, equivalently, the Jaynes- 
Cummings model), which describe the dipole- field inter- 
action, may be simplified by adiabatically eliminating the 
emitter polarization, so that the emission kinetics in such 
a system is adequately described by the traditional laser 
rate equations [7[: 

^= P(i)-7||A^(i) -PlwNmt) (1) 
^= -TMt)+(3l\\N{t) +/37||A^(t)s(t) (2) 

where P{t) is the incident pumping rate; s{t) is the num- 
ber of photons in the lasing (reference) mode; l/Tc and 
1/711 are respectively the lifetime of the reference mode 
(also termed cavity lifetime) and the dipole lifetime; and 
/3 is the fraction of spontaneous emission funnelled into 
the reference mode. 

These are continuous-variable differential equations 
that describe the evolution of the mean photon and dipole 
populations, but ignore the "granularity" , that is the dis- 
crete nature of these populations. In any particular real- 
ization of the lasing process, the photon and dipole popu- 
lations can assume only integer values, with discrete tran- 
sitions between the different population configurations. 
The physical processes behind these transitions are the 
excitation of a dipole through "pumping" , the emission of 
a photon by a dipole inside or outside the reference mode, 
or the escape of a photon from the reference mode to the 
outside. The discrete nature of the photon and dipole 
numbers can be taken into account rigorously (and un- 
der the same assumptions of rapid loss of phase memory) 
through the master equation approach, which considers 
the probability of each realization of the discrete photon 
and dipole populations, as well as the transitions between 
the different population configurations 0. It can there- 
fore account for the full statistics of the emission process. 
Within the master equation approach the discrete tran- 
sitions between different population configurations can 
be regarded as probabilistic partitions of the photon and 
dipole numbers produced by the processes of pumping, 
emission, and escape, governed by the statistical param- 
eters of the master equation. 

When the master equation is applied to the calculation 
of the mean photon and dipole numbers ((s) and (A^)), 
it yields the traditional laser rate equations ^ and 
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under the assumption that the fluctuations in the dipole 
andphoton numbers are uncorrected {{Ns) = {N){s)) 
0, S], which is a good approximation when the num- 
ber of photons and dipoles is large (so that their relative 
fluctuations are small) or, more precisely, in the "ther- 
modynamic limit", that is when 1//3 oo. 

Calculation of the laser photon statistics by use of the 
master equation requires calculating at each time-step 
a matrix of {Nmax x Smax) interdependent probabilities 
(where the subscript "max" indicates a cut-off value cho- 
sen large enough so as to encompass all population values 
with a significant probability), a formidable task when 
the dimensionality of the matrix is larger than a few tens. 
Alternatively, the master equation may be solved through 
stochastic simulation methods [1, [l^ , whereby the sys- 
tem evolves from one discrete state to another through a 
jump event in each time-step, thus describing a trajectory 
in state-space. The average of a large number of trajec- 
tories converges to the solution of the master equation. 
This greatly simplifies the calculation, as each trajectory 
is calculated independently of the others, and conver- 
gence of the average is reasonably rapid. Beyond the 
mathematical convenience it offers in the calculation of 
the master equation, the calculation of individual trajec- 
tories permits a direct visualization of the time-evolution 
of the system parameters and thus provides physical in- 
sight into the ongoing processes. 

The laser rate equations provide a guideline in the cal- 
culation of the stochastic state-space trajectories, since 
these equations deal with the mean values that may be 
calculated through the master equation. Thus, at each 
step, the mean of the stochastic evolution of a given ini- 
tial state must follow the laser rate equations. 

To this end, the rate equations dD) and © can be 
solved numerically by combining an iterative procedure 
over a time-step dt, small compared with l/Fc and 1/7||, 
and a quantum jump approach to ensure that N and s 
are integers. The use of a small time-step ensures that 
this iterative procedure follows well the fine structure of 
the time-evolution of the system, while at the same time 
it preserves the correlations between the fluctuations in 
the photon and dipole populations. Thus, re-writing the 
laser rate equations, the discrete evolution of N and s 
follows at each iteration step: 



N{t + dt)= N (t) + P{t)dt- 
s{t + dt) = s(t) - Tcs(t)dt +{1 + 



+ l3s{t))N{t)dt (3) 
sit))pj^\Nit)dt (4) 



This presentation of the laser rate equations permits us 
to interpret each increment or decrement of N and s 
as resulting from a random drawing from an appropri- 
ate distribution, such that its average is the correspond- 
ing expression in these equations. Thus, denoting by 
PD[^] a random number out of a Poisson distribution 
with mean value /x, and by BD[n,p] a random number 
drawn from a binomial distribution of n trials, each of 
which yields success with probability p, we can obtain 
the number of additional dipoles excited by the pumping 
process during time-step dt as PD[P{t)dt], while Nd{t), 



the number of dipoles having decayed between t and 
t + dt hy spontaneous or stimulated emission, is given by 
Ndit)=BD[N{t),-/\\ (1 + l3s{t))dt]. Each photon emitted 
by these dipoles has a probability {1 + s{t)) (3 / {1 + (3s{t)) 
to be tunneled in the cavity mode, and thus the num- 
ber of emitted photons entering the cavity mode is 
BD[Nd{t), (1 -I- s\t))l3/{l + I3s{t))]. Finally, the number 
of photons that escape to the outside due to cavity losses, 
is drawn from the binomial distribution BD[s{t)^Tcdt\. 

This procedure accounts for the evolution of the pho- 
ton and dipole populations as integer numbers. It thus 
takes account of the "discretization noise" without any 
need to introduce it in the continuous-variable laser 
rate equations externally. Traditionally, noise (including 
quantum fluctuations) is re-introduced in the laser equa- 
tions (H)) and ^ by linearizing around a (large) mean 
value of the populations and adding a Langevin driving 
force with Gaussian statistics [l^. However, this ap- 
proach is valid only in the case of large populations [15j , 
where the relative fluctuations are small and obey Gaus- 
sian statistics thanks to the central limit theorem. 



III. POPULATION OSCILLATIONS IN 
NANOLASERS 

To gain an insight into the effects of discretization 
noise, we have calculated state-space trajectories at dif- 
ferent pumping powers, by solving Equations ([3]) and ^ 
numerically through the above procedure, using the val- 
ues of ^ = 0.25, 1/7|| = 400 ps and l/F^ = 10 ps (cor- 
responding to a quality factor of 13000 at A = 1.5/im), 
which are typical for quantum dot nanolasers [iHSj]. 

Figure 1 (left) displays the time-dependence of the 
number of intra-cavity photons and excited dipoles at a 
pumping power of A.bPthi where Pth = Fc//? is the nom- 
inal threshold pumping power. At this pumping power, 
the mean number of dipoles has reached its pinning value 
of Fc//37|| ~ 160, yet the photon and dipole populations 
display well-defined oscillations of the same frequency, 
as can be seen in the corresponding power spectra (Fig. 
1 right). These oscillations are synchronized (but in 
quadrature), in a manner analogous to the population 
cycles driven by demographic stochasticity observed in 
finite predator-prey ecosystems [Tsj . 

The synchronization of the photon and dipole popula- 
tion cycles as a function of pumping power can be studied 
through the covariance of these two quantities, as can be 
seen in Fig. [5] (left) (black continuous line), where the 
covariance curve for a low-/? laser (red continuous line) is 
also given for comparison. Below threshold, the photon 
and dipole populations fluctuate randomly and indepen- 
dently of each other, as the presence of less than one 
photon (on average) in the cavity produces no change 
in the dipole emission dynamics and spontaneous emis- 
sion outside the cavity mode destroys any correlation be- 
tween these two populations. At higher powers, as stim- 
ulated emission grows, a larger fraction of the emission 



4 




1,0x10- 



2,0x10-' 
Time (s) 



3,0x10-' 



4,0x10-" 




10" 10" 10'" 10" 

Frequency (Hz) 



10'^ 



FIG. 1: (Left) Evolution of the number of photons s{t) (blue line) and the number of excited dipoles N{t) (red line) as a 
function of time for a nanolaser with /3 — 0.25, 1/7|| = 400 ps and l/Fc = 10 ps for a pumping power of P = 4.5Pth, with 
Pth = Tc//?, the nominal threshold pumping power. (Right) The corresponding power spectra. 
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FIG. 2: (Left) Covariance of the number of photons s{t) and the number of excited dipoles N{t) as a function of the pumping 
power, for /3 — 0.25 (black line) and /3 — 0.001 (red line). Dotted curves are calculated through the relaxation oscillation theory. 
Curves are rescaled by the pinning value of the dipoles Fc//37||. (Right) Comparison of the photon population cycle frequency 
for P = 0.25 (continuous line) with the relaxation oscillation frequency (dotted line), as a function of pumping power. 



goes into the reference mode, so that any decrease in the 
dipole population is matched by a synchronous increase 
of the intra-cavity photon population, thus producing a 
strong negative correlation. The difference between high- 
and low-/3 lasers is essentially the speed with which the 
rescaled covariance reaches its limiting value of -1. At a 
given pumping power (normalized with respect to thresh- 
old), a high-/? laser has a smaller number of intra-cavity 
photons, so that its stimulated emission rate is smaller 
implying that a smaller fraction of the emission goes into 
the reference mode. As a consequence, the rescaled co- 
variance drops more gradually to its limiting value of -1. 

The power dependance of the covariance (Fig. [2] left) 
and of the central frequency of the population oscillations 
(Fig. fright) are quite similar to those calculated through 
the traditional linearized ocillation relaxation theory [l9[ 
(Fig. [2] dotted lines), indicating that the photon and 
dipole population cycles in a nanolaser can be under- 
stood also as arising from relaxation oscillations driven 



by the noise due to the granularity of photon emission. 
This mechanism is similar to that of the origin of the 
excess linewidth in semiconductor lasers UM and of the 
excess noise in class-B (i.e. Tc > 7||) lasers |18i|, however 
with an important difference. In traditional low-/3 lasers, 
the (spontaneous and stimulated) emission events are dis- 
tributed essentially uniformly in time, corresponding to 
a white noise spectrum that gets filtered by the transfer 
function of the system to produce a damped oscillation. 
On the other hand, in high-/3 lasers, because of the small 
mean value of the intra-cavity photons, emission events 
are more frequent at the top of the photon population 
cycle, and thus constitute a periodic driving force that 
produces a resonant amplification sustaining the photon 
and dipole population cycles. 

The photon population cycles should be directly ob- 
servable in the second-order photon correlation function 
(7^^-'(r), which is measurable experimentally through a 
Hanbury Brown and Twiss setup. As can be seen on Fig. 
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FIG. 3: Second-order photon correlation function g^'^\r) of 
a nanolaser (/3 = 0.25) undergoing photon population cycles, 
for P = Q.hPth (red line) and P = 4.5Pth (black line). The 
dashed horizontal line corresponds to g^^^r) = 1, the value 
in continuous-wave lasers. 



mode. However, the modulation was weaker at a pump- 
ing current of 2.72 mA (11 % above threshold and ~ 1500 
photons in the lasing mode) and was below instrumental 
sensitivity at currents of 3.00 mA (23 % above threshold 
and ~ 2300 photons in the lasing mode). It was at- 
tributed by the authors to relaxation oscillations excited 
by (generic) noise in the laser. 

We argue that in that experiment the oscillations are 
driven by discretization noise, as our model can fit the ex- 
perimental second order autocorrelation function (7^^)(t), 
with no adjustable parameters other than the cavity life- 
time l/Fc = 3.4 ps and the spontaneous emission lifetime 
l/7|l = 4 ns extracted from [17|. Figure |4] shows the ex- 
perimental second order autocorrelation function ^'^■'(t) 
(black points) for a pump current of I = 2.6 mA (cor- 
responding to a pumping power 6.5 % above threshold 
and ~ 1000 photons in the cavity) and the fit using the 
discretization noise model is represented as a continuous 
red line. 



[3] (red line), below threshold g'^^^r) displays the usual 
behavior for chaotic light. That is, when the two de- 
tected photons are in temporal coincidence it has the 
value S^^HO) = 2, while when the two photons are sepa- 
rated by a time interval r the value of g^^^ (r) drops to 1 
exponentially e~^'''^ , where Fc is the cavity energy decay 
rate. This is the familiar Siegert relation, 



g(2)(r) = l+ .gW(^) 



(5) 



where g^'^'^ir) is the first-order correlation function and 
corresponds to the Fourier tranform of the optical spec- 
trum. On the other hand, above threshold, the presence 
of photon population cycles in nanolasers causes the value 
of g^'^\Q) to be larger than 1, in contrast to traditional 
lasers which display g'^^^ir) = 1 for all time intervals r. 
At finite time-intervals, ^'•^■'(t) tends towards the value 
of 1, in a damped oscillation fashion dropping periodi- 
cally below 1 for a few cycles (see Fig. [3] black line). This 
feature arises from the presence of population cycles in 
nanolasers (even under constant pumping) and is a devi- 
ation from the Siegert relation (Eq. (5)), which does not 
pemit g^^^T) values smaller than 1. 



IV. OBSERVATION OF POPULATION 
OSCILLATIONS 

In a recent publication [l7| , observation of a fast peri- 
odic modulation in the photon correlation g^^^r) was re- 
ported for a single-mode vertical-cavity surface-emitting 
semiconductor laser displaying a spontaneous emission 
coupling factor of /? « 10~^ (as can be seen on Fig. 2 of 
Ref. [l3|)- The modulation was clearly observable at a 
pumping current of 2.47 mA, corresponding to a power 
of 1 % above threshold and 150 photons in the lasing 
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FIG. 4: Second order autocorrelation function q'-^^(r) as 
function of time. Black dots: experimental data from [ITI ]. 
Red continuous line: simulation through the discretization 
noise model. 



V. POPULATION OSCILLATIONS AND NOISE 

For small-/? lasers (in which the number of photons 
at threshold is large), or when the nanolaser response 
is averaged over times much longer than the relaxation 
oscillation period (and thus corresponds to an ensemble- 
averaged response), the population cycles arising from 
granularity appear as fluctuations around the (mean) 
steady-state populations. These fluctuations are often 
characterized by the Fano factor [18;], 



F : 



s 



s{g('H0)-l) + l 



(6) 



where is the variance in the photon number of mean 
value s. By analogy to phase transitions, it is often ex- 
pected that fluctuations should be maximum at threshold 
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FIG. 5: (Left) Power dependence of the Fano factor for /? — 0.25 (black line) and /3 = 0.001 (red line). (Right) Power 
dependence of the normalized root-mean-square noise for /3 — 0.25 (black line) and /3 — 0.001 (red line). The dotted line (of 
logarithmic slope -1) is a guide to the eye. 



0,123], as can be seen in Fig.[5]for a low-/? laser. However, 
for class-B lasers (F^ > 7|| ) with large /3, the Fano factor 
maximum occurs at powers significantly higher than the 
stimulated emission threshold (see Fig. [5] left), as shown 
by van Druten et at. [Tsj . While this discrepancy be- 
tween the laser threshold and the fluctuation maximum 
has been interpreted as an indication that such lasers 
produce partially chaotic light "21^ and become coherent 
only at pumping powers much higher than that of the 
Fano maximum, our results indicate that the increase of 
the Fano factor not due to the presence of a chaotic com- 
ponent in the laser output but arises from the large am- 
plitude of the photon population cycles that are excited 
and sustained by discretization noise. 

An alternative designation of the magnitude of the fluc- 
tuations, useful in engineering, is the normalized root- 
mean-square noise (nRMS), which is related to g^'^\Q) 
by 



nRMS = 



= Wg(2)(0)-l + 



(7) 



As can be seen on Fig. [5] (right), for traditional low-/? 
lasers and above threshold, the nRMS noise drops rapidly 
to very small values. As the signal-to-noise ratio is essen- 
tially the inverse of nRMS, this feature makes the output 
of low-/3 lasers an ideal carrier for high-fidelity informa- 
tion transmission. For high-/? nanolasers, on the other 
hand, the nRMS noise remains quite high even at pump- 
ing powers much above threshold, both because of the 
presence of population cycles (which cause a large value 
for ^'■^•'(O)) and the relatively small value of s. As the 
signal-to-noise ratio is of the order of 1 even high above 
threshold, this implies that such lasers are ill-adapted to 



the transmission of information, unless particular noise- 
reducing protocols are implemented taking into account 
the population cycles. 

VI. CONCLUSION 



In conclusion, the output of a continuously-pumped 
high-/? nanolaser cannot obey the traditional (continu- 
ous, real-number) laser equations, because of the small 
numbers of intra-cavity photons and dipoles involved, 
which bring forth the discrete nature of these small 
populations. Such nanolasers should display population 
cycles both for the intra-cavity photons and the excited 
dipoles because of the intrinsic stochasticity of the 
photon and dipole granularity, analogous to what is 
observed in small ecosystems. These cycles can be 
understood also as resulting from relaxation oscillations 
driven synchronously by the sudden and discrete jumps 
in photon and dipole numbers, due to the elementary 
emission processes. The photon population cycles are 
observable as a periodic modulation of the second-order 
correlation function of the laser, as was recently reported 
by Takemura et al. 



17| . The presence of population 



cycles in high-/? nanolasers may appear as excess noise 
which may hinder information processing and transmis- 
sion with such lasers. 
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